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Abstract. In this paper we study existence and nonexistence of nonnegative distributional solutions 
for a class of semilinear fractional elliptic equations involving the Hardy potential. 
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Introduction 

Let B be a ball of R N , N > 3, centered at 0. Let s £ (0, 1), p > 1 and 7 > 0. In this 
paper, we study existence and nonexistence of nonnegative functions u £ C].r\L^ oc {B) 
satisfying 

(0.1) (-A) s u--f\x\- 2s u = u p in 5, 

where 

Equality (0.1) is understood in the sense of distributions. The distribution (— A) s u 
is defined as 

((-A) s u,ip) = [ u(-A) s ipdx V<peC™(B). 
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Here the fractional Laplacian (— A) s is defined via the Fourier transform as 

(0.2) (-A)V(x) = / \C\ 2s me t( - x dC, 

(27t)t Jrn 

where 

£(C) = ^)(0 = -^-jr [ e~<-Mx)dx 

is the Fourier transform of ip. The nonlocal structure of the fractional Laplacian 
(— A) s can be seen in its representation in the real space: 



(0.3) (-A)V(z) = C S , N [ 



<f{x) - <p(y) 

\ X - y \N+2s 



dy, 



for some positive constant C Sj at. For the equivalence between (0.3) and (0.2), we 
refer the reader to [32]. 

Problem (0.1) is related to the Hardy inequality which were proved by Herbst in 
[28] (see also [45]): 



(0.4) 70 / \x\- 2s u 2 dx < [ \(\ 2s u 2 d( Vu G C, 

Jr n Jr n 

where 

p2 ( N+2s \ 

70 = 2-^^^. 



(N-2) 



The constant 70 is optimal and converges to the classical Hardy constant , 
when s — >■ 1. Here T is the usual gamma function. We should mention that (0.4) is 
a particular case of the Stein and Weiss inequality, see [41] . 

A great deal of work is currently been devoted to the study of the fractional 
Laplacian as it appears in many fields such as probability theory, physics and math- 
ematical finance. We refer the reader to papers [10], [39], [40], [24], [4] (and the 
references there in) for a nice expository. A good reference for the potential theory 
of (— A) s can be found in the book of Landkof [32]. The operator (— A) s — 7o|a;| _2s 
appears in the problem of stability of relativistic matter in magnetic fields. One 
can see [38] where a lower bound and a Gagliardo-Nirenberg-type inequality were 
proved. 
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The problem of existence and nonexistence of (0.1), for s = 1, was studied by 
Brezis-Dupaigne-Tesei in [6] where the authors showed that for /3 G [0, ^y^] , 1 < 
P < jv-2-2/3 ' ^ ne problem 

_An - (^ (iV ~ 2)2 - ^ |x|~ 2 = u p in V'(B) 

has a positive solution u G L P (B) and does not have any nonnegative and nontrivial 
supersolution u 6 Lf oc (B \ {0}) when /3 G [0, ^j-^) an d P > jv-2-2g • Some related 
results and problems are in [5], [6], [15], [17], [18], [23], [43], [21], [22], [20], [7], [8]. 
Our results in this paper extends the one of Brezis-Dupaigne-Tesei in [6] to the case 
s G (0, 1). Before stating them, we fix the following notation: for a G [0, N ^ 2s ), we 
put 

r (N±2s±2a\ r ( N+2s-2a \ 
= 2 2s L \ 4 ) L \ 4 ) 
< a y ^ N-2s-2a ^ y ^ N-2s+2a ^ " 

The mapping a H> ^ a is monotone decreasing and ^ a — > when a — > N ^ ls ■ We 
should mention that this constant appears in the perturbation of the 'ground-state' 
| a; | 2 for the operator (— A) s — 7o|x|~ 2s . Indeed, letting i9 a = \x\~ hQ we have 

(-A)^ a - 7Q |x|- 2 ^ a = inl^\{0}, 

see Lemma 3.1 in Section 3. 

Our existence result is the following 

Theorem 0.1 Let a G [0, ^=^] and 1 < p < j^zffEff • There exits a function 
u G C\ n L P {B) satisfying u > in B and 

(-A) s u--f a \x\- 2s u = u p inV'(B). 

As what concerns nonexistence, we have obtained: 

Theorem 0.2 Let a G [0, ^=^) and u G C\ n £f oc (£ \ {0}) stic/i that u > and 

(-A) s u- la \x\- 2s u>u p inV'(B\{0}). 
IfP> N-ll-la ' thenu = 0mB. 
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We observe that if a = we have p + 1 = j^% s ■ the critical Hardy-Littlewood- 
Sobolev exponent and that ^-ll-ia ~^ +°° as a ~^ ^T^- 

The proof of Theorem 0.2 relies on weak comparison principles recently used by 
the author in [23]. However substantial difficulties have to be overcome due to 
the nonlocal structure of the fractional Laplacian. Nonexistence result of nonlinear 
elliptic problems using comparison principles have been obtained in [1], [2] [37], [35], 
[31], [36], [30] and the references therein. 

For the existence result, in the supercritical case N-2s-2a ^ P — N-2s > we nave an 
explicit solution constructed via In the subcritical case, p+ 1 < > we use d 
standard variational arguments thanks to the following improved fractional Hardy 
inequality: 

Theorem 0.3 Let 2 > q > max ^1, j^^j ■ Then there exists a constant C > such 
that for all u G C™(B), 

(0.6) C\\u\\ 2 w r, q{B) < 70 / \x\- 2s u 2 dx- [ \C\ 2s u 2 d(, 

v ' JB JR N 

where r = — ^ . 

This result, which might be of self interest, is proved in Appendix 5. 
The proof of all the results presented above are manly based on a Dirichlet-to Neu- 
mann operator B s for which B s u = (— A) s u in D'(E) for any Lipschitz bounded 
open set E of R N and u is the null extension outside £ of a function u belonging 
to some Sobolev space. To be more precise let us first recall the result of Caffarelli 
and Silvestre. We recall that 

H S (R N ) = {u : R N -> R : (1 + \(\ s )u € L 2 (R N )} . 

Put 

= {( t ,x) : t > 0, x e R N }. 

Given w £ H S (R N ), minimization procedure yields the existence of a unique function 

If 



T-L(w) G ff 1 (]R^ r+1 ;t 1 2s ) being the harmonic extension of w over the half space 



(0.7) 



&v(t l - 2s VH{w)) = in R% +1 , 
T-L(w) = w on 
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In [12], Caffarelli and Silvestre proved that (— A) s w is given by the Dirichlet-to- 
Neumann operator lim^o t l ~ 2s 9n @^ : 



lim i-2 S d?M. = Ks (-A) s w in R N , 
t^o dt v ' 



for some constant k s > 0. In addition 

(0.8) I t 1 - 2s \VU(w)\ 2 dxdt = k s I |C| 2s ^C- 

Jr^ +1 Jrn 

We want to provide similar arguments in bounded open sets. We define the 
Hilbert space S> S ' 2 {R N ) which is the completion of C£°(R N ) with respect to the 
norm: 

«->• / \(\ 2s \v\ 2 d(- 
Jr n 

Let E be a bounded open set in R N with Lipschitz boundary. We introduce the 
Hilbert space 

J^ S (E) := {u G H S {E) : u G @ S ' 2 (R N )}, 

where 

_ u in E 
u = < 

in R N \ E. 
The space J^f(E) is endowed with the natural norm 

(o.9) \\ u \\ 2 = [ ici 2 -|5| 2 dc=/ |ci 2 W)l 2 *;. 

v ; Jr^ Jrjv 
Note that, since E is bounded, by (0.4) there exists a constant C(E) > such that 

(0.10) C(E)\\u\\ H s {RN) < ||«||jr '(E) < I|«I|jt«(R") Vn G JT S (£). 

From this we deduce that 

(0.11) ^o(E) = {u£ H S (E) : u G H S (R N )}. 

Hence, see for instance [[27], Theorem 1.4.2.2 ], the space C^°(E) is dense in r J%f s (E). 
By (0.11) for any u G Jt? s (E), we can consider its harmonic extension %{v) as in 
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(0.7). We define the Dirichlet-to-Neumann operator B s : J% S (E) ->• J^~ S (E) given 
by 



t-+o dt ' 

where J%? ~ S (E) is the dual of Jt? s (E). This operator turns out to be linear and it 
is an isometry, 

\\Bsu\\jp-°(E) = \\u\\jr s (E) € 3#q{E) 

by (0.8). Moreover 

B s u=(-A) s u mV'(E) Vu€Jf '(E). 

In particular a solution u G Jt? s (E) to the problem 

£> s ii = / in E 1 

yields a solution v G ^ S ' 2 (R^) to the problem 

(-A) s t> = / in £7, 

u = in R N \ E 

and conversely. We refer to the next section for more details. 

In order to get, say, qualitative informations on the solution to the problem 

B s u=(-A) s u = f, 

it is, in general, more convenient to work with the (mixed) problem 

div^-^VWCu)) = in R£ +1 , 

New difficulties arise here because of the (possible) degeneracy of the equation (0.12). 
However the weight t 1-2s falls into the Muckenhoupt class of weights thus regular- 
ity results, Harnack inequalities are available (see [19]) and this is enough for our 
purpose in this paper. 

An interesting characterization of J#q(E), see [27], is that Jt? s (E) is the interpola- 
tion space (H$(E),L 2 (E)) S;2 : 

' H S {E) a € (0,1/2), 
(0.13) ^o(E) = l H$ (E) 8 = 1/2, 

k Hq(e) se (1/2,1), 



(0.12) 
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where 



i 



H^E) = \uem{E): J e ^dx < oo 



endowed with the natural norm, with d{x) = &ist(x,dE). 

Remark 0.4 Let E be a smooth bounded domain ofR N . Recently a pseudo differ- 
ential operator A s of order 2s was introduced by Cabre and Tan [11] for s = 1/2 
(see [14] for every s / 1/2J in the following way: for any u G J^ S (E) 

oo 

k=l 

where fik is the zero Dirichlet eigenvalues of —A with corresponding orthonormal 
eigenf unctions ipk and uu = f E u(fkdx is the component of u in the L 2 (E) basis 
Wk}- 

Using (0.13), it was shown in [11] and [14] that 

(0.14) J^ s (E) = ^ueL 2 (E) : M |K| 2 < ooj . 

The operator A s corresponds to the Dirichlet-to-Neumann operator given by the har- 
monic extension over the cylinder E x (0, oo). Indeed, let H\ S (E x (0, oo)) be 
the set of measurable functions w : E x (0,oo) — > E with w G H 1 ^ x {r\,r2)), 
< r\ < r2 < oo and w = on dE x (0, oo) such that the following norm 

Him (£x(o,oo))= / t^lVwfdxdt < oo. 

In [11] and [14], the authors showed that for any g G J^~ S {E) there exists a unique 
solution u G Jt? s (E) to 



(0.15) 



A s u = g in E, 
u = on dE. 



In addition u is the trace of w G H\ S (E x (0,oo)) which is the unique solution to 
(0.16) 



div(t 1_2s Vu;) = in E x (0, oo), 
w = on dE 1 x (0, oo), 
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where kn, s (kn,s = 1 f or s = 1/2 j is a constant depending only on N and s. 
Moreover, it holds that, with the norm in (0.14), 

(°- 17 ) IMP = K ^IMI^ s (£x(0,oc))- 

We can compare the operator A s with the operator B s . For simplicity, we consider 
the case s = 1/2. Assume that g G C^°(E) and u is a solution to (0.15). Take w 
its extension over the cylinder. Consider H(u) which is the harmonic extension of 
u in IR^ +1 given by (0.12). Clearly 



H(u)>w inR^ +1 . 

It follows from Hopf lemma that 

dw &H (u) ^ 

-m > —^ L mE - 

Hence 

A 1 / 2 u > Bi/ 2 u in E. 

In particular the operator A s yields (up to a multiplicative constant) subsolution 
to the fractional Laplacian (— A) s . This is the reason why the use of B s is more 
convenient in this paper. 

We give here the plan of the paper: 

• Section 1: Notations and Preliminaries. 

• Subsection 1.1: Dirichlet-to- Neumann operator. 

• Section 2: Comparison and maximum principles. 

• Section 3: Nonexistence of positive supersolutions. 

• Section 4: Existence of positive solutions. 

• Appendix 5, Subsection 5.1: Remainder term for the fractional Hardy inequal- 
ity. 
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1 Notations and Preliminaries 



Let u G L 2 (M. N ), we will consider its Fourier transform 

£(0=^(0(0:=—^ [ e-<- x u(x)dx. 

For s > 0, the Sobolev space H S (R N ) is denned as 

H S (R N ) = {ue L 2 {R N ) : \(\ s u G L 2 ^)} 

with norm 

IMItf^R^) = II^IIl 2 (IRJV) + |||C| S: "|| L 2( R iV). 
We also have by Parseval identity 

\ u{x) - u(y)\ 2 
\x - y\ 



I«IIh.(r") - ll-"lli 2( iR^) + / / ' dxd y- 



Let E be a bounded domain in with Lipschitz boundary. For q > 1 , we introduce 
the space W s,q (E) defined as the space of measurable functions u such that the 
following norm is finite 

wi • wi • f f ~ u (y)\% a 



\W«(E) — U"ULi(E) ^ J E J E \ X -y\N+g 

We define Wq' q (E) to be the closure of C%°(E) with respect to the norm || • || y/ a , q ^)- 
As a notation convention, we put H S (E) = W 1 ' 2 (E) and H$(E) = Wq' 2 (E) which 
are Hilbert spaces. 

It is well known that if u € H l (E) then u, its null extension outside E, is in 
// 1 (IR Ar ) and = ||ii||iji(M Ar ) ■ This is not in general true for functions in 

H S (E) (s = 1/2 for instance). We shall define a space of functions in which we 
recover this defect by imposing integrability of null extensions. 

The Hardy inequality (0.4) suggests the definition of the Hilbert space & s ' 2 
which is the completion of C^°(M N ) with respect to the norm: 



(1.1) I \(\ 2s \v\ 



'■d(. 



9 



u = < 



As it will be apparently clear in the remaining of the paper, we introduce the Hilbert 
space 

(1.2) ^o( E ) := {u G H S {E) : u G $> S ' 2 (R N )} , 
where we put here (and hereafter) 

u m.E 
in R N \ E. 
The space Jif s (E) is endowed with the norm 

(1.3) \\u\\\s {E) = I \C\ 2s &\ 2 dC= [ \t\ 2s \T(u)\ 2 dC 

Note that, since E is bounded, by (0.4) there exists a constant C{E) > such that 

(1.4) C{E)\\u\\ H s^ N) < \\u\\ M ^ {E) < \\u\\ H .^ Vn G J^o( E )- 
Therefore 

je o s (E) ={ue H S (E) : u G H S (R N )}. 
See for instance [[27], Theorem 1.4.2.2 ], the space C™{E) is dense in J% S (E). 

Notations : For G an open set of R N , we use the standard notations for weighted 
Lebesgue spaces: L p (G;a(x)) = {u : G — > R : f G u p a(x)dx < oo} . B N (0,r) is a 
ball in R N centered at with radius r > and S^ 1 = dB N (0, 1). R* +1 = {(t, x) : 
t>0,xe R N }. B^ +1 (0,r) = R^ +1 n B N+1 (0,r) and = R^ +1 n S N . If there 
is no confusion, we will put B N = B N (0,1) and B^ +1 = #f +1 (0, 1). The space 
W${B% +1 (p,r)-,a(x)) = {u G ^(5^ +1 (0,r);o(i)) : u = on S^}. 



1.1 Dirichlet-to-Neumann operator 

It is well known that the space of Schwartz functions <S contains C£°(R N ) and that 
J 7 is a bijection from S into itself. In particular (-A)V G C°°(R N ) n L°°(R N ) for 
every if G CfflR^). In fact we have for any if G C 2 {R N ), see [39], 



(-A)V(x)| < C jtff^ VxgrJV - 



This motivates the following: 
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Definition 1.1 Let G be an open subset of W N . Given u G C\, the distribution 
(-A) s u G V(G) is defined as 

((-A) s u,ip)=[ u(-A)"<pdx V^GC C °°(G). 
Jr n 

Some recent results conserning s-superhamonic functions in the sense of distributions 
as above are in [39]. 

Consider the Poisson kernel of := {(t,x) : t > 0, x G R N } 

(1.5) P(t, x) = PNtS t 2a 



\ x 2 +t 2|(JV+2a)/2' 



where pn, s is a normalization constant, see [10] for an explicit value. Let 
we can define 

u(t, x) = P(t, •) * u = p N>s t 2s [ U{V) N+2s dy V(t, x) G R?+i_ 

JRN (| y _ x |2 + t 2)^ 

It turns out that 

div(^- 2s Vu) = 

Therefore u is smooth in R^ +1 . Moreover if u is regular at some point xo then 

limn(t,xo) — > u(xo). 
By an argument of [12], we have that 

(1.6) -limt 1 - 2 ^^) = Ks (-A) s u(x ), 
where the constant k s is explicitly computed in [10]: 

(1.7) *. r(1 - s) 



2 2s (l -s)T(s) 



For any w G H S (R N ), we denote by T~L(w) its unique harmonic extension over R^ +1 . 
Namely (see for instance [12], [10]) U{w) G i? 1 (R^ +1 ; t 1_2s ) and 

div(^- 2s VHH) =0 in R^ +1 , 

(1.8) \ U{w) = w on R^, 

_ tl - 2s dH(w) = Kg (_ A y w on M iv. 
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In particular if w G C 2 (R N ) then %{w) = P(t, •) * w. In addition one can check (see 
[12] ), using integration by parts and the Parseval identity, that 

(1.9) / t 1 - 2s \V'H(w)\ 2 dxdt = k s f \(\ 2s wd( = Ks [ \{-&) s/2 w\ 2 dx. 

Jr* +1 JR N Jr n 

Therefore from the definition of the space r J%? s (E), we have 

(1.10) K a \\v\\%,s (E) = [ t l - 2s \VH(v)\ 2 dxdt VvGJ^iE), 

o JrN+i 

where as usual v is the null extension of v outside E. 



We now introduce a Dirichlet-to-Neumann operator B s defined on J#q(E) which 
can be seen as a local version of (1.8). 

Proposition 1.2 Let E be a bounded open set with Lipschitz boundary. Denote by 
JT- S (E) the dual ofJ^ s (E). Then the mapping B s : J^ S (E) -»• J^- S (E) given by 

(B s v, <p)je-'{E)^e '(E) = - K 7* J } i f tl ' 2s ^r'^ dx Vv > V G -*%( E ) 
is a linear isometry. In addition for any v G J4? S (E) we have 
(1.11) B s v = (-A) s v inV'(E). 

Proof. By definition for any v G J^ S (E), v G H S (M. N ) thus the operator B s is well 
defined and linear. Consider %{v) which satisfies (1.8). Then integration by parts 
yields for every ip G J#q(E) 

[ t^VHiv) ■ VUmdxdt = [ hint 1 

j^N+i j RN t ^o dt 

This, (1.10) and Holder inequality imply that 

(B s v,ip) j#>-s(E),je£(E) — IMI^(£)IMIjr s (,E) 

while 

(B s v, «)jf- S (E) ; jr s (E) = ll v lljr s (£)- 

Hence 

ll^s«||jr- s (E) = IMI,^ s (£)- 



,l-2s^( v ) ~j 

* - -ipdx. 
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On the other hand, for any tp G C%°(E), we have by integration by parts 



[ t l - 2s VU{v)-VH($)dxdt = [ 

- f 



limt 1 



dt 



-tpdx 



t-5>0 



dt 



This means that 



= k s / v(—A) s ipdx. 
Jr n 

B s v = {-A) s v mV'(E). 



□ 



We turn to the characterization of the space (E) . As suggested with the fact 
that H{v) G H 1 (R r ^ +1 ;t 1 - 2s ) for every v G J% S (E), we have the converse: 

Proposition 1.3 Let E be a bounded open set with Lipschitz boundary. Define 



Hl T (E-t l - 2s ) = |u> G H 1 



+ 



1 ;t 1 ' 2s ) : w 



= on 



\e} 



and 



Hfc(E) = {u G V S ' 2 (R N ) : u = in R N \ E) 
We have the following equalities: 



(1.12) J4f s (E) = { 

In particular 



u\ : u G 
\e 



Hfc(E)} = [w\ e : w G Hl T (E;t 1 ~ 2 ')} 
(-A) s u = B s u inV'(E), Vu€^ s T (£), 



where u = u 



E 



Proof. The first equality in (1.12) is immediate by definition. The second equality 
is a consequence of the trace embedding theorem. Indeed, take w G H^ T [E]t l ~ 2s ). 
Then the null extension of w outside E is nothing but w which belongs to H S (R N ) 

E 

and in addition ||w||# s (£;) < IMI_fp(RJV)- D 
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Summarizing, we state the following 



Proposition 1.4 Pick g G 3f S (E). Let v G 3%q(E) (given by the Lax-Miligram 
theorem) be the unique solution to 



Let w £ H 1 



1. +l-2s 



B s v = g in E. 
) solve the mixed problem 

div(t 1-2s Viu) = inR£ +1 , 
w = onR N \E, 



-t 



l-2sdw 
dt 



K s g on E. 



Then v = w in E; for any ip G 3f s (E) 



(BsV,ip)jr- s{E) ^ s{E) = J |CI T(v)T(<p) 



= f {-A) s ' 2 v{-A) s l 2 ^dx 
Jr n 

= kT 1 [ t^VwVUmdxdt 

= k- 1 [ i 1_2a VW(v) • Vn(v)dxdt 
Jr n+1 



and thus 



I t'- 2s \v 



w\ 2 dxdt. 



K s\\ v \\,je Q a {E) 

We can extend the above in unbounded domains: 

Remark 1.5 Here we consider E any open subset of~R N with Lipschitz boundary. 
Define 



(1.13) 



3T S (E) := {u G H S (E) : u G H 1 ^)} , 



where as usual u stands for the null extension of u outside E. We have that C£°(E) 

is dense in 3t s (E), see [27]. 

By similar arguments, we have that the operator 



B 8 (V) = -K-H 



dt 



+ v 
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is a linear isometry form J^ S (E) -»• (Jf s (E))' , where (Jf s {E))' is the dual of 
je s (E). 

2 Comparison and maximum principles 

Unless otherwise stated, E is a bounded Lipschitz open set of W N . We have the 
following technical result which will be useful in the sequel. 

Lemma 2.1 Let E n be a sequence of Lipschitz open sets such that E n CC E n+ \ and 
U^i-EW = E. Let g n G L 2 (E) such that g n -> g in L 2 (E). Consider v n G J$? Q s (E n ) 
solution to 

B s v n = g n in E n . 
If v G J4? S (E) is the unique solution to 

B s v = g in E 

then Un — ^ v in L 2 {E). 

Proof. Observe that U{n n ) G H^E-t 1 - 28 ) thus by Proposition 1.3 G J% S (E). 
In addition we have by Hardy and Holder inequality 

II^IIj^OE) = lkn||jT s (E) < C(E)\\g n \\ L 2( E y 

Therefore Un is bounded. By assumption it converges weakly to v in Jif s (E) and 
strongly in L 2 (E) because C™(E) is dense in J^ S (E). □ 

The following maximum principle can be found in [[14] Lemma 2.4] or in [19]. 

Lemma 2.2 Let E be a bounded Lipschitz domain ofR N . Let v G J^ S (E), v > 
such that 

B s v > in E. 
If v then for any compact set K C E 

ess inf v > 0. 

K 
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Lemma 2.3 Let g G L 2 (E), g > and let w G Lj oc (R+ +1 ), such that 

/ t l - 2s \Vw\ 2 dxdt < oo 
Jr^ +1 

and 

(2.1) / t^^Vw -V<pdxdt + c wtpdxdt>K s gcf>dx 
Jr^ +1 Je Je 

for every <j) G T (E; t 1_2s ), where c > 0. Assume that w > on R-^ \ £\ T/ien 
w > in rJ TT . 

Proof. Test (2.1) with max(— w,0) G H^E;^- 23 ). □ 

Lemma 2.4 Let c > and Zei n G C\, u > and g G L 2 (E) such that 

(2.2) (-A) s u + cu> 5 inV'(E). 
Let v G Mq{E) solves 

(2.3) fi s v + ci> = # in £\ 
T/ien 

u > v in E. 
Proof. Recall that (2.2) is equivalent to 

(2.4) f u(-A) s ipdx> [ g<pdx-c[ mpdx V<p G C^(E), <p > 0. 
Jr n Je Je 

Denote by p n the standard mollifier (which is symmetric: p n (—x) = p n (x)) and put 

U n = Pn* U. 

Claim: for any ip G 



(2.5) / (-A) s u n cp= [ n(-A) s ( Pn *^). 

JR N JR N 

It is easy to check using Fubini's theorem and the symmetry of p n that 

(2.6) / (-A) s u n (pdx = / u n (-A) s (pdx = up n * (-A) s ipdx. 
Jr n Jr n Jr n 
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Now we notice that, in M. N , 

Pn * (-A)V = HHPn * (-A)V)) = HK\ 2S (HPn)Hv))) = {-^) S (Pn * <P). 

Using this in (2.6), we get (2.5) as claimed. 

Let E n := {x G E : dist(x,&E) > 1/n}. We deduce from (2.4) and (2.5) that for 
all ip e Cg°(E n ) and ip > 

/ (-A) s u n (pdx = / u(-A) s (p n * ip)dx > / g(p n * ip)dx - c u(p n *(p)dx 
Jr n Jr n Je Je 

= (Pn* g) i fdx - C (p n * u)<pd,X. 

Je Je 

We conclude that 

(2.7) (-A) s u n (x) + cu n (x) > p n * g(x) =: g n (x) for every x G E n . 

Pick X G C™(R N ) such that x = 1 on E n . We let w n = H(xu n ) € fl' 1 (R+ +1 ;t 1-2s ) 
be the harmonic extension of x^n: 



(2.8) 



It turns out that 



l-2s 



div^-^Vwn) = R£ +1 , 



duo 

(2.9) -t 1 ' 23 -*^ 1 + cti; n = K s (-A) s u n + cu n > K s # n on £„. 

Let v n £ Jif s (E n ) be the solution to 

B s v n + cw n = g n in £ n . 



Then it follows that w n > %{v n ) in R^ +1 by Lemma 2.3. In particular u n > v n in 
£" n . By Lemma 2.1, — >■ t> in L 2 (E) and the proof is complete. □ 

We recall the definition of the s-capacity of a compact set A C E: 
(2.10) C S (A) = inf {II^H^^b) : <P > 1 in a neighborhood of A}. 

<f>£C£°(E) 

Note that if C s (^4) = then \A\ = by Poincare inequality (see (1.4)). We have the 
following comparison result modulo small sets. 
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Lemma 2.5 Let A be a compact subset of E with C S (A) 
g G L 2 (E) such that 

(-A) s u + cu> g inV'(E\A). 

Let v G J^ S {E) solve 

B s v + cv = g in E. 

Then u> v in E. 

Proof. Let A £ be a smooth open e-neighborhood of A compactly contained in E. 
Define D £ = {x G D : dist(x, d(D \ A £ )) > e}. It is clear that 

(-A) s u + cu> g mV'(D e ). 

Consider v £ G J^ S (D £ ) solving 

B s v £ + cv £ = g in D £ . 

By Lemma 2.4 we have u > v £ in D £ . The same argument as in the proof of Lemma 
2.1 yields v £ G M^{E) for every s G (0,1) and it is bounded. Hence it converges 
weakly to some function w in J4? Q S (E) and strongly in L 2 (E). In particular u > w. 
Moreover for any (p G C£°(E \ A), we can choose e > so small that suppcp is 
contained in D £ thus taking the limit as e — > 0, we get 

( w ^)^(e) + c j E WL P = J E g(pdx G c Ti E \ A )- 

From this equality, to conclude the proof (that is v = w), it suffices to show that 
C™{E \ A) is dense in C™(E) with the J^ s (£)-norm because w G J%?(E). 
Since C S {A) = 0, there exists a sequence ip n G C^°(E) such that ^„ > 1 in a 
neighborhood of ^4 and in addition 

(2- 11 ) IIXn||^( E ) < IIV>n||^(£) -> 0, 

where \ n = min(^ n) l)- Now take any <f> G C^°(E) and note that (1 — Xn)4> £ 
C c °°(£ \ A) and moreover (1 - Xn)</> -> ^ in ^ s (£) by (2.11). This concludes the 
proof. □ 



= 0. Let u G C\ and 
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We shall define a new space which is more convenient when dealing with the 
Hardy potential. Namely, we assume that there exists b G Lj oc (E) and a constant 
C > such that 

(2.12) \W\\\s{e) ~ j Hx)^ 2 dx >cj <p 2 dx G C™{E). 

Definition 2.6 Let b G L\ oc (E) so that (2.12) holds . The Hilbert space Jf * b (E) is 
the completion ofC^(E) with respect to the scalar product 

Note that the Lax-Miligram theorem implies that for any / G L 2 (E), there exits a 
unique solution to the problem 



(2.13) 



B s v — b(x)v = f in E, 
in the sense that for all <fi G J^ s b (E) 



( v ' $)^{E) ~ J b{x)v<j)dx = J f<j)dx. 



Remark 2.7 Let e > 0. Put d £ {x) = b(x)(l - e). Then J^ S {E) = ^ ds {E) by 
Propositon 1.3. This holds true because if v G ffl§ d {E) then by (2.12) we have 
veH s (R N ). By similar argument J^ S {E) = J%* b {E) if b G L°°(E). 

Lemma 2.8 Let A be a compact subset of E with C S (A) = 0. Let b G Lj oc (E) such 
that (2.12) holds. Suppose that u G C\ with u,b>0 and f G L 2 (E) such that 

(2.14) (-A) s u- b{x)u > f inV'(E\A). 

Let v G Ji? s b (E) be the unique solution to 

B s v — b(x)v = f in E. 

Then 

u > v in E. 
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Proof. For every r > 0, we define E(r) = rE and we put 

A(r) = inf MjgCgW) > 

<p£Cg°(E(r)) \\(p\\ L 2(E(r)) 

Notice that A(r) = r _s A(l) for every r > 0. Without loss of generality, we may 
assume that A(l) > 1. Fix k G N, k > 1. We rescale by setting 

= 7^TT7 u (^2fe x ) ' 6jfc(ar) = 6(\/2fear) , f k (x) = f(V2kx), A k = -^=A 
V>k) s y/2k 



Note that (-A)V* = (2fc) s (-A)V for every <p G C™(E) with ^(x) = p(V2fca;). 
Thus u k satisfies 



(-A) s u k - b k {x)u k > f k in V'{E(1/V2k) \ A k ) . 
For any fixed k, we let vq G J^ s (E , (l/v / 2fc)) be the solution to 

B s vo = f k in E(l/V2k). 
Then wo < in E(l/y/2k) by Lemma 2.5. We truncate the weight b k (x) by letting 

b k {x) = mm(b k (x),k 2s ). 
For ii g N, we let v n G J% s (E(l/V2k)) be the solution to 

B s v n = h(x)v n ^i + / fc in E(l/V2k). 
We use once again Lemma 2.5 to infer that 

vo < ■ ■ ■ < v n -i < v n < u k in E(l/V2k) Vn G N. 
In addition by Holder inequality, 

\\ v 4j^(E(i/V2k)) ^ A(l/V2fc)"" / ' fc " L2( ' E(1/V ^ )) 

and by induction 

n k 2is 

\\ V n\\^f(E{l/VZk)) ^ Wfk\\ L z(E(i/V2k)) Zj \(1/V2k) 2i+1 



- ((2*)^))=' 

oo ^ 

- H/fell^^l/v^)) Z 2^ < 00 " 
i=0 
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We conclude that v n — 1 v k in J4? S (E (1/V2k)) asn->oo. In particular v k solves 
(2.15) B s v k = b k (x)v k + f k inE(l/V2k) 

and v k < u k in E(l/y/2k) (because J4? S (E(1/ \f2k)) is compactly embedded in 
L 2 (E{l/\fZk))). We scale-back by letting 

v k (x) = 2k v k (-^x 



'Ik 

so that v k < u in £, v k G ^ S (E) and moreover by (2.15) for all ip G C%°(E) 

(2.16) (w fe , p)^ a {E) ~~ / min(6(x), k 2s )v k (pdx = / /(^dx. 
Hence 

ll« fc H*» fc (i5) ^ ll^ll^(B) ~ J^mm{b(x),k 2s } \v k \ 2 dx 
= j^fv k dx<C\\v k \\^s b(E) 

by Holder inequality and by (2.12), where the constant C depends on / but not on 
k. Therefore the sequence v k is bounded in Jtf s b (E). We conclude that there exists 
v G J4? s b (E) such that, for a subsequence, v k — ^ u in J4? s b (E). Now by (2.16), we 
have 

^ k ^)je YE) + I (Kx)-^n{b(x),k 2s })v k i p= I ftp. 

J E J E 

Since 

(6(x) - min{6(a;), A; 2s }) € fc p < (6(x) - min{6(x), fc 2s }) u|p| < 2&(x)u| ¥ >| G L l {E), 
the dominated convergence theorem implies that 

(2.17) (v,p)^ b(E) = f f<p forany^GC-^). 

J E 

We conclude that v = v by uniqueness. By (2.17), we have 

kn2 je a b {E) = J^ b (E) ~ ( v ^ k )jff s b (E) + {v,V- ^ k )^l h {E) 

= (^^ k )^ b (E) ~ fv k + (V,V~ V k )jg> o s b{E y 
E 



\v — V 



21 



Now notice that 

je> b {E) < ll« fc H3r '(E) " / ram{b(x),k 2s }\v k \ 2 dx = ! fv k . 

J E J E 

From the two above inequalities, we deduce that 

C j \v - v k \ 2 dx < (v,v - v k )j^s b ^ -)• 

by (2.12). Hence v k — > v pointwise and thus v < u in E. □ 

Remark 2.9 The same result as in Lemma 2.8 holds if we assumed the coercivity 
that there exist constants C, c > such that for all ip G C£° (E) 

(2.18) \\<p\\ %>s^ + c / tp 2 dx — / b(x)(p 2 dx>C / ip 2 dx. 

J E J E J E 

One works with 



^\M\s {E{r)) +c\\p>\\ 2 L2{E{r)) 

A c (r) = ml - — 

V dC™{E(r)) \m\tf(E(r)) 

instead of \(r) and notes that A c (r) = r~ s A(l) + c. 

We close this section with the following useful lemma and its immediate consequence. 
Its counterpart, for s = 1, is in [22]. 

Lemma 2.10 Let E be a bounded Lipschitz domain of "U. N . Let A be a compact 
subset of E with C S (A) = 0. Let u G C\, b 6 Lj oc (E) and u, b > 0. Assume that 

(2.19) (-A) s u>b(x)u inV'(E\A). 
Then 

(2.20) / t l - 2s \VU{p)\ 2 dxdt = k s \\v\\%s (e) > k s [ b(x)<p 2 dx \/ip G C^°(E). 

Proof. Put gk(x) := min(b(x)u, k) > for integers k > 1. Let Vk G Jt? s (E) be the 
solution to 

B s Vk = 9k in E. 
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By Lemma 2.2, we have ^ € L^ C (E) and by the standard maximum principle 
T~L{vk) > 0. Moreover by Lemma 2.5, we have 



(2.21) u > v k > mE. 



Let ip e C™(E). Put V k = U{v k ) and V k £ = V k + e, for e > 0. Set V = W so that 
V k £ tp 2 G #q T (.E; ^~ 2s ). Simple computations show that 

|w%)| 2 = | v fc £ wf + w fc £ • v(y fc e v 2 ) = |v?v?/f + w fc • v(v£ v 2 )- 

Thus using integration by parts we have 

„2 



Je 



9k ~ | 32 < 



Take the limit as e — > to get 

I t l - 2s \VU{<p)\ 2 dxdt > k s I 9 —^ 2 dx 

Jr^ +1 JEVk 

by Fatou's lemma. By (2.21), we infer that 

[ t l - 2s \Vn(^)\ 2 dxdt >k s [ 9 —ip 2 dx. 
Jr^ +1 Je u 

Again by Fatou's lemma, inequality (2.20) follows immediately by taking k — > +oo. 
□ 



Theorem 2.11 Let E be a bounded Lipschitz domain o/IR^ with £ E, N > 3. 
Then there is no nonnegative and nontrivial u G C] satisfying 

(-A) s u > ~/\x\- 2s u inV'(E\{0}), 

with 7 > 7o- 

Proof. Note that C s ({0}) = provided N > 2s. If such u exits then u > in E 
by the maximum principle thus Lemma 2.10 contradicts the sharpness of the Hardy 
constant 70. □ 
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3 Nonexistence of positive supersolutions 



We start with the following 



Lemma 3.1 For every a G (— ^ — s, ^ — s), put # a (x) = \x\ 2 +a . Then 
(-A) s $ a = la \x\- 2s $ a inR N \{0}, 

where 

(3.1) Ja 



p ^ N+2s+2a ^ p ^ jV+2s-2q -j 
" p ^ N-2s-2a ^ p ^ iV-2s+2a ^ 



For a > 0, the function a 1— >■ 7 Q is continuous and decreasing. 
There exists a positive function T a G C' 3 \ {0}^ suc/i £/ta£ 

div^ 1 - 25 ^) = in R* +1 
(3.2) J T a = $ a on dR^ +1 \ {0} 

~ tl ~ 2sQ ¥ = « S (- A )^a = ^7aM" 2S tfa On C^ +1 \ {0}. 
Moreover if a > t/ien |VT a | G L 2 (5^ +1 (0, i?); t 1 " 28 ) /or every R > 0. 

Proof. Note that i9 a G The Fourier transform of radial functions (see 
Theorem 4.1]) yields 

1-iV f°° 1 iV-1 

J\ya)\P)=P 2 / (r/3) 2 J jy-2 (rp)v n r 2 rfr, 







where J iv-2 is the Bessel function. Then we have 

2 



/-oo 

F(^ a )(p) = P — — / M s+a j^MdM 

Jo 2 



where 



'0 

m Q/ 9 2 



p f iV+2 J s+2a N \ 

m a = 2 s+a -^ * ^ 



p ^ jV-2s-2a^ - 

Now we notice that 



(-A) s tf Q = F(F((-A) s tf Q )) = F(p 2s F($ a )(p)) = m Q F(p^+^ a ) = la r- 2s $ 
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T a (t,x) = < 



For the proof of the fact that the map a H> 7 a is continuous and decreasing, we refer 
to [16]. 
We define 

' P(t,-) *# a (x) V(^)eRf +1 

Mx) Vx G dR* +1 \{0}, 

where P is the Poisson kernel defined in Section 1.1. Clearly T a is positive. We 
have that 

Fit 

-t 1 ' 2 ^ = Ks (-Ay# a inR"\{0}. 

Hence we get (3.2). 

From the regularity theory of [19], we deduce that T a G \ {0}) for some 

P > 0. In addition T a G H 1 (fi x t 2 ); t 1 " 25 ) for every fi CC 1^ \ {0} and 
< ti < t 2 < oo. 

Finally if a > we have that $ Q is the trace of the function = +a 

which satisfies 

-div(t^ 2s y a ) = ( (iV ~ 2 ' )2 - a 2 ) t 1 " 25 !^- 2 ^ > in 

where z = (x, i). It turns out that T a < V a . From this, we deduce that t 1_2s |z|~ 2 T Q , G 
L 2 (B^ +1 (0,i2)) and also |x|~ 2s T a G Lf oc (R N ) because a > 0. We let <p be a cut-off 
function such that <p = for \z\ < e, <p = 1 for 2e < |z| < i?, <^ = for \z\ > 2R and 
\Vtp\ < Ce" 1 for e < \z\ < 2e. We use (p 2 Y a as a test function in (3.2) to get 



-2st2 „2 



/ ^- 2s VT Q V(^ 2 T a ) = / |x|- 2s T 2 

Integrating by parts and using Young's inequality, for some constant c > 0, we have 
cf ^- 2 V 2 |VT Q | 2 < [ \x\- 2s Tl^+[ t 1 " 2 ^ 2 ^ 2 . 

JB% +1 (p,R)\B% +1 (p,e) JB N (0,2R) Jr^ +1 

Therefore 



cf t'- 2s \VT a \ 2 < f \x\- 2s T 2 a +f t l - 2s \{t,x)\~ 2 T 2 a . 

JB^ +1 (0,R)\B* +1 (0,e) JB n (0,2R) Je<\(t,x)\<2e 

Fatou's lemma yields |VT Q | G L 2 (B^ +1 (0, R); t 1 - 28 ). □ 
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The comparison result obtained in Lemma 2.8 allows us to derive the following 
estimate when the potential b(x) is the Hardy one. 

Lemma 3.2 Let N > 3, a G [0, (TV - 2s)/2) and p > 1. Suppose that u & £}, Pi 
Lf oc (B N (0,2) \ {0}), u > sucft f/ioi 

(3.3) (-A) s n-7 Q |x|- 2s u > n p mD'(/(0,2)\{0}). 

T/ien there exists a constant C > suc/j £/taf 
(3.4) 

\\<P\\le'(B»(o 2)) - 7« / |x|- 2 V 2 <te > C [ <p 2 dx V<p G C?(B N (0, 2)). 

' V ' " JB N (0,2) JB n (0,2) 

Moreover there exists a constant C > suc/j i/iai 

(3.5) « > u Q > C'lxl 2 ^^ in 5^(0,1), 

where v a G Jt? s b (B N (0, 2)) fioitfc 6(x) = 7 a |x|- 2s j is i/te solution to 

B s v a - ~i a \x\- 2s v a = min(n p , 1) in B N (0, 2). 

Proof. Inequality (3.4) is trivial for a > so we consider only the case a = 0. 
Since C s ({0}) = provided TV > 2s, by Lemma 2.5 and Lemma 2.2, we have u > 
in 5^(0,2) and 

M = ess inf u > 0. 

5^(0,1) 

Hence it G satisfies 

(-A) s u - j a \x\- 2s u > M p ~ l u in V'(B N (0, 1) \ {0}). 

By Lemma 2.10 we obtain (3.4) thanks to the scale invariance of the integrals on 
the left hand side. 

We put b(x) = 7 Q |x|- 2s . By (3.4), for a > 0, we can let v a G J% s b (B N (0,2)) be the 
solution to 

B s v a - ~i a \x\- 2s v a = min(u p , 1) in B N (0, 2). 

Then by Lemma 2.8, we have u > v a in B N (0, 2). We first consider the case a > 0. 
Then 7 a < 70 and thus v a G (5^(0, 2)). By the regularity result of [19], we 
get that V a = H(v^) is continuous in B^ +1 (0, 1) \ {0} and V a > by Lemma 2.3. 
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Consider $ a and its harmonic extension T Q given by Lemma 3.1. The maximum 
principal (see Lemma 2.2) implies that we can set 

min v a 

S N 

(3.6) C = > 0. 

max I a 



Put w = C'T a — V a . We have weakly 
(3.7) 



dw^-^Vw) = in B^ +1 (0, 1) 
w < onSf, 

01 



< W < 



2S W -W a \x\- 2s w < on 5^(0,1). 



Then u; + := max(u>,0) G 5 (i?+ +1 (0, 1); t 1 2s ) and therefore by integration by 
parts 

f t l ~ 2s \Vw + \ 2 dxdt - k s7q f \x\- 2s (w + fdx < 0. 
Jr^ 1 Jb n (0A) 



'(0,1) 
In particular 

\\ w+ \\%s(B"(o,i)) ~ 7a / |xr 2s (w; + ) 2 da; < 0. 

o v \ , )i Jb n (0,1) 

Hence w + = by Hardy's inequality. Hence v a > C'-da in B N (0, 1) that is (3.5) for 
a > 0. 

For the case a = 0, we put ot n = 1/n and we notice that the sequence v an 6 
^(£^(0,2)) solution to the problem 

^a„-7a„N" 2 X„ = / 111^(0,2) 

is monotone increasing to vq because the mapping a t— > ^y a is decreasing. Therefore, 
taking into account (3.6), we readily get (3.5). □ 

Proof of Theorem 0.2 

Lemma 3.3 Let E be a bounded Lipschitz domain of M. N , N > 3. Suppose that 
€ E and a G [0, (N - 2s)/2). Let u G C\ n £f oc (£ \ {0}) sucfc i/iai 

(3.8) {-A) s u-~f a \x\- 2s u > u p inV(E\{0}), 
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with 



la = 2 



p2 I N+2s+2a ^ 
2s L \ 4 



p2 ( N-2s-2a ^ " 

Proof. Assume that u / 0. It follows from Lemma 3.2 that there exist r,C r > 
such that 

(3.9) u{x) > v a (x) > C r |x|^ +a Vx G 5^(0,7-) C E, 
where v a G Jt? Q s b (E) (with 6(x) = 7 Q |a;| _2s ) is the solution to 

•Bs^q — 7a\x\~ 2s v a = min(n p , 1) in E. 
On the other hand Lemma 2.10 yields, for all ip G C^°(E), 

(3.10) \\f\\%»(E) ~ 7a / kr 2 V 2 ^> / vT^dx. 

Je Je 

We first consider the case a > 0. If r is small, by (3.9) we have, for < a 1 < a, 

(-A) s u- la ,\x\- 2s u> {- lal+la + cr 1 )\x\~ 2s u inP'^COUO}). 

By Lemma 2.8 and using the same arguments as in Lemma 3.2 we get, provided 
a' /■ q, 

(3.11) u{x) > C' r \x\^ K+a ' VxG B N (0,r/2), 

for some constant C' r > 0. Using the estimate (3.11) in (3.10) we get 



2 

Jt? »(B N (0,r/2)) 7a 



\x\~ 2s ^dx > (C' r f- 1 / |x|(^+ a ')(p- Vdz, 

B N (0,r/2) JB N (0,r/2) 



for all <p G C™(B N (0,r/2)). Since p > ^§5^, we have 

-5 := (^^- + c/) (p - 1) + 2s < 0. 
Hence for every p G (0,r/2) 

^ W 0,p)) > (7a + KY^P- 5 ) [ \x\- 2s <p 2 dx y V G C?(B N (0,p)). 
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This contradicts the sharpness of the Hardy constant thanks to scale invariance of 
the inequality. 

Finally, for the case a = we note that (3.10) implies, by density, that 
\\ v a\\\s(E)-l<* j \x\- 2s v 2 a dx> I v p a +l dx. 

J E J E 

This also leads to a contradiction because v a G Jt? Q s b (E) while by (3.9) 
/ v p+1 dx > C f \xf-^)^dx 

JE JB N {Q,r) 

> C f [ t~ l dtda = +oo, 
Js N - 1 Jo 

for some constant C > 0. □ 

4 Existence of positive solutions 

The proof will be separated into several cases. We put 

E a (u) := |MI|^(B) - 7a J \x\~ 2s u 2 dx Vu G C™(B), 

where B is a ball in M. N centered at with N > 3. 

Case 1: a € (0, (N - 2s)/2] and 1 < p < (N + s)/(N - 2s). 

Thanks to the Hardy-Littlewood-Sobolev inequality and Hardy's inequality we have 

/ , \2/(p+l) 

(4.1) E a (u) > U b uP+1 ) Vu G C T(B). 

Thanks to the compact embedding of J^ S (B) into L P+1 (B), we can minimize E a 
over the set 

(4.2) LeJ? Q s (B) : J B (u + r +1 = l\. 

Let u G J^q(B) be the minimizer. Put it 1 * 1 = max(±u, 0). By Proposition 1.3, 
belongs to J4? S (B). We check rapidly that E a (u + ) < E a (u). Observe that 



JR f f +1 JK N+1 
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because H(u ± ) and ^{(u)^ 1 have the same trace on R N while %{u ± ) has minimal 
Dirichlet energy. Now using this and Hardy's inequality we have 

K s E a (u + ) = K s \\u + \\ 2 % ,s {B) - K s J a / \x\~ 2s (u + ) 2 dx 

J B J B 

<? II +II 2 / I l~ 2s _l_ II ~ll 2 

— ^s||^ ll^ s (_B) K s^fa J | X | ti Ota; -+- K s 1 1 U \\j^(b) 

< f t l - 2s \V(U(u)+)\ 2 dxdt-K sla I br 2 Vdx 

Jr1 +1 Jb 

+ / t 1_2s |V(H(«)~)| 2 dx^ 

= [ t 1 -*\VH(u)\ 2 dxdt-K s 'y 0l [ \x\' 2s u 2 dx 

Jb 

= K s E a (u). 

Thus u = u + . Therefore there exists a Lagrange multiplier A > such that 

B s u — j a \x\~ 2s u = \u p mB. 

i 

Hence \p- 1 u is a solution of problem (0.1). 
Case 2: a = and 1< p < (N + 2s) /(N - 2s). 
Lemma 5.4 yields for every q G ^2, max ^1, j^^j^j 

E (u)>\\u\\ 2 w r, q{B) uGC^(B), 

with r = - J. Therefore Jf s >b (B) is compactly embedded into L p+l (B), with 
& = 7oM _2s - Hence we can minimize Eq over the set 

(4.3) : J B (u+y +1 = iy 

We have to check again that £'o(n + ) < Eq(u). But this can be done by density and 
using similar arguments as above. We skip the details. We get a positive minimizer 
u = u + of £b in the set (4.3). We conclude that \~u is a solution to (0.1) for some 
Lagrange multiplier A > 0. 
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Case 3: a G (0, (JV-2a)/2) and (7V+2s)/(iV-2s) < p < (N + 2s - 2a)/(N -2s- 2a). 
Consider dp = r~ 1-/3 given by Lemma 3.1 which satisfies 

(-Ay#p = 1 p\x\- 2s #p mR N \{0}. 

-2s 

We look for a solution of the form w = [irp- 1 with a constant \i > to be determined 

2s — N | a —2s 

in a minute. Assume that we can take /3 > such that r 2 +p = rP -i then 

(-A) s w = ~ip\x\~ 2s w + w p - w p ~ 1 w 

= j a \x\- 2s w + w p + (79 - 7 a - /i p_1 )|x| _2s u;. 

Since /3 i-)- 7/3 is decreasing, we can choose = 73 — 7a > provided a> ft. But 
note that a > /3 as soon asp < (iV+2,s-2a)/(iV-2s-2a) andp > (N+2s)/(N-2s) 
implies fi > 0. In conclusion we have, in \ {0}, 

(-A) s u> - <y a \x\~ 2s w = w p 

and to G n L p (£) is a solution to (0.1). □ 



5 Appendix 

5.1 Remainder term for the fractional Hardy inequality 

Let E be a bounded open set of 1^, iV > 3, with £ E. The following (local) 
Hardy inequality is a consequence of (0.4) 

(5.1) 70 / u 2 \x\~ 2s dx < \\u\\%,s (E) Vu£C™(E). 

J E 

In addition the constant 70 is optimal. Our objective, in this section, is to improve 
inequality (0.4) in bounded domains of WL N . 

Many deal of work has been done in improving the classical Hardy inequality start- 
ing from the work of Brezis- Vazquez [9]. We also quote [3], [44], [26] for related 
improvements. 

We shall prove a Vazquez-Zuazua-type (see [44]) improvement for the fractional 
Hardy inequality (0.4). That is for 2 > q > max ^1, 2~^j , there exists a constant 
C(E) > such that for all u G C£°(E), 

C{E)\\u\\ 2 w r, q{E) < ||n||^ s(E) -70 / |x| a_ Vcte, 
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where a = 1 — 2s and r = — - . The proof requires several preliminary lemmata. 

Consider the function T defined in Lemma 3.1 satisfying 

div^-^To) = inM^ +1 , 
(•)-2) <(T = |x|T „„ R N {()}. 

.^1-2,^0 =Ks7o | x |-2, To onM JV\ {0} _ 

We have seen in the proof of Lemma 3.1 that 

(5.3) |T (*)| < C\z\( 2s ~ N V 2 Vz = (x,t) G B^ +1 . 
By scale invariance, we have that 

(5.4) T {z) = R {N - 2s)/2 T {Rz) Vi?>0. 
This implies the estimate 

(5.5) \? (z)\ > C\z\^-^ 2 VzGBf 1 
and also 

(5.6) |VT (s)| < Clzl^-W-i Vz€B» +1 . 

We now prove the following result which were proved in [16] when s = 1/2. 

Lemma 5.1 For every q G (1,2) there exists a constant C > such that for all 
if G C?(B N+1 ) 

/ \ 2/q 

(5.7) C\ f &\V(p\ q dz) < f t a \V^\ 2 dz - k s1o f Ixl^^dx, 

\Jb^ +1 J Jb^ +1 Jb" 

where a = 1 — 2s. 

Proof. Let ip G C%°(B N+1 \ {0}) and put tp = Simple computations yield 

|V^| 2 = |T VVf + vr -v(ToV 2 ). 

Integration by parts and using (5.2) leads to 

f t a \V<p\ 2 dz - k s7o [ |a;| a - Vdc > f t a Tl\V^\ 2 dz. 
Jb^ +1 Jb" Jb^ +1 
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By (5.5) and using polar coordinates z = ra = Mjff, we get 

(5.8) f t a T 2 \V^\ 2 dz >C f [ (a^rlVi/^fdadr, 

Jb1 +1 Jo Js^ 

where o\ is the component of a in the t direction. We wish to show that there exists 

a constant C > such that 

(5.9) 

/ l \ q / 2 
I -If [ (<7i)°r|W| 2 d<Tdr | >cf t^\Vip\ q dz \/<p e C™(B N+1 ). 

We have 

f t^\V<p\ q dz= [ t^\V^T + ^VT \ q dz < C [ t 3 ^ (\V^T \ q + \i;VT \ q )dz. 
Jb^ +1 Jb^ +1 Jb^ +1 

Put 

I 1= f if |V^T | 9 cfe, 
Jb» +1 

Using (5.3) and Holder inequality, we have 

h < C I (a,)^ 2 C ' r N+qi - l ' N ^ 2 \V^\ q drda 
Js* Jo 

(5.10) = C f 1 r K +q (l-N)/2- q /2 [ (<Tl )W2 r 5/2| V ^|«d(7dr 

Jo Jsy 

1 / \ q l 2 

< C r N{2 ~ q V 2 (J Ja^rlV^daj dr 



< C {^J N {a{) a r\V^\ 2 d^j dr 

< C ^ J Ja^rlV^dadr^j 



< CI. 
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On the other hand we have using (5.6) 

h< C ! (a^ 2 f ' r N ^' 2 -^\ q drda 
Js? Jo 



< C f {ai) qa/2 f 1 r N(2-q)/2-q/2+q 

Js" Jo 



dip 



dr 



Q 

drdo 



(5.11) < C f 1 r N+^-N)/2- q /2 f ( ai )W2 r? /2| VV) | ?d(jdr 

Jo Js^ 
< CI, 

where in the second inequality we have used the one dimensional Hardy inequality 
/J f q dr < c Jq r~ q \f'\ q dr and observing that (5.10) is just (5.11). The lemma follows 
because C™(B N+1 \ {0}) is dense in C™(B N+1 ) with respect to the H 1 (B N+l ;t a )- 
norm when N > 3, see [29]. □ 

The Lion's interpolation inequality, [[33] Paragraph 5], shows that for a / and 
— 7, < % < h there exits a constant C > such that 

q z q 

(5.12) 

c \\ v \\wr 9m ^< f t^\Vv\ q dxdt+ [ t^\v\ q dxdt, Vw G C c 0O (R JV+1 ), 
w { > Jr^ +1 Jr^ 1 

where r = ^=2 - -. 

2 q 

We first prove a generalized weighted Poincare trace inequality. The proof is stan- 
dard. We recall that the space W 1 J(B+ +1 ;t 9 °/ 2 ) was defined in Section 1. 

Lemma 5.2 Suppose that q > 1 andr := 2 -5p — | > 0. T/ien the following inequality 

holds 

(5.13) 

/ i ^ a / 2 M 9 (fe <C f i t qa / 2 \Vu\ q dz + C I \u\ q dz Vu G W^J(5^ +1 ;t W2 )- 

Proof. Inequality (5.13) is well known for a = 0. So we restrict ourself to the case 
a + 0. 

Assume by contradiction that (5.13) does not hold. Then there exits a sequence 
u n G Wl' q {B^ +l ;t qa ' 2 ) such that 

(5.14) / t qa l 2 \Vu n \ q dz+ [ \u n \ q dz = o(l) 

Jb n+1 Jb n 
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and 



f t qa l 2 \u n \ q dz > VneN. 
Jb" +1 



Up to normalization, we may assume that J b n+i t qa / 2 \u n \ q dz = 1. But then (5.14) 
implies that u n is bounded in W l ' q {B™ +1 ; t qa ' 2 ) thus u n -± u in W^-B^ 1 ; t«°/2) 
and u„ -> « in L q (B^ +1 ;t qa / 2 ) (see [25]) so that 



(5.15) f t qa ' 2 \ 



u\ q dz = 1. 



It follows from (5.12) and the compact embedding of W^ q {B N ) into L q {B N ) , that 
ti n — ► u in L q (B N ). From (5.14) we get it ^ = and also Vu = 0. It turns out 
that u = in -B+ +1 , a contradiction with (5.15). □ 

By an argument of partition of unity, we have the following 

Lemma 5.3 Let 2 > q > max ^1, 2~^) ■ There exist some constants C, c > swc/t 
that for all (p€C^(R N+1 ) 

(5.16) C||v9||^ T , 9fBJV s < / t a \Vip\ 2 dz - K s ~f \x\ a ~ 1 ip 2 dx + c ip 2 dx, 
( ' Jr^ +1 Vr^ Jr^ 

where a = 1 — 2s and r = — 7- . 

A q 

Proof. We put 

J(v)= f t a \Vv\ 2 dz- k s7o f \x\ a ~ 1 v 2 dx. 

JR^ 1 Jrn 



Let x G C c °°(5 Ar+1 ), < x < 1 in and such that x = 1 on B N+1 (0, 1/2). Let 

iJV- 

+ 



rj G // 1 (i?^ r+1 ; t a ) be the minimum of the problem 



inf 

Then 



y^ +i t a \Vu\ 2 dz : <z- X e 

div(t a V7?) = in B^ +1 
r\ = x on B N 
r] = onSf. 
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It turns out that < rj < 1 in B+ +1 . In addition, thanks to [13], lim t ^ t a ^ £ 
L^ C (B N ). Given ip G C%°(R N+1 ), simple computations based on integration by- 
parts lead to 

I t a \V(pr])\ 2 dz < [ t a \Vip\ 2 dz + c [ ip 2 dx, 

where c > depends only on rj. On the other hand we have 

/ {x^iriipfdx = I \x\ a - 1 i P 2 dx+ I (1 - 77 2 )|x| a -V 2 
Jb n Jr n Jr n 



< / \x\ a ~ 1 ip 2 dx + c / ip 2 dx. 
Jr n Jr n 



Therefore we obtain 



Applying Lemma 5.1, we infer that 

(\ Vi 
[ t%\V(ri<p)\ q dz) <J{f)+c[ f 2 dx. 
Jr1 +1 J ! Jrn 

By Lemma 5.2 and Holder inequality (1 < q < 2) 

C\\w\\ 2 w r, q(K+ ,^ tqa/2) < J(<p) + cj^ V 2 dx. 
Using Lions' interpolation inequality (5.12) with a / 0, we obtain 

CWvvWwT.qm*) <J((p) + c ip 2 dx. 

JR N 

If a = 0, it is well know that W n ' <? (R+ +1 ) embeds continuously into W l - l / q > q (M. N ). 
Recalling that r\ = x = 1 on B N (0, 1/2), the lemma follows by scaling. □ 

Taking advantages to the singular nature of the Hardy potential and the scale 
invariance, we prove the main result in this section: 

Lemma 5.4 Let 2 > q > max ^1, ■ Then there exists a constant Cq > such 
that for all ueC™{B N ), 

(5.17) Coll 

u llw-g ' q (B N ) — \\ u \\'jf s (B N ) ~ To J \ x \ a ^u 2 dx, 
where a = 1 — 2s and r = 2 ^ L — \ . 
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Proof. Let u G C™(B N ) and we define U = H(u) = H(u). Then U G t a ) 
thus by Lemma 5.3 and a density argument we get 



C\\U\\wr q(R N^ < I t a \VU\ 2 dz- k s7o [ \x\ a ~ l U 2 dx + c / U 2 dx. 



\W t >i(B n ) 

Since u = U on 5^, it follows that 



c \\ u \\w-r, q(B N\ < \\u\\%>s (b n) ~ 7o / M a 1 u 2 dx + c u 2 dx. 

Let r G (0, 1) we derive from the above that for every u G C^°(B N (0, r)) 
(5.18) 

CIMIw'-.ifB^fO.r)) ^ H U ll^(B^(0,r)) "TO / |x| 0_ Vda; + C / U 2 cfe, 

V V " v v " JB N (0,r) JB N {0,r) 

with c, C > independent on r. This holds because |MIj^ s (-6-^(0,7-)) = ll' li ll l ^ ) s (s iV (o,i)) 
as long as u G C™{B N (0,r)) and r < 1. 

It is clear from (5.18) that wee need only to show that there exists a constant 
C 2 > such that for every u G C™{B N ) 

(5.19) C2 / u 2 dx < \\u\\ 2 yiPam N\ — In 1 \x\ a ~ 1 u 2 dx. 

J B N ^0(« ) J BN 

Before proceeding, we recall that the mapping a i-> j a (defined in Lemma 3.1) is 
decreasing. We will use this fact and the estimates in Lemma 3.2 to conclude the 
proof. Pick 

/ N-2s 
a G 0, — 



and let r > be so small that 

(5.20) ( 7o - la )r a ~ l - c> 0. 



As we did in Section 4, by (5.18), we can define the space Jf? s b (B (0,r)) with 
b(x) = 7o|x| a ~ 1 — a Letting 2 < p + 1 < j^ N 2s , we can choose q (close to 2) so that 
WQ' q (B N (0, r)) is compactly embedded into L P+1 (B N (0, r)). Then minimization 
procedure implies that there exits a nonnegative and nontrivial u r G Jt? s b (B N (0, r)) 
solution to 

B s u r + cu r - 7o|x| a ~ 1 n r = Cu v r in B N (0, r). 
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By density 

(5.21) (-A) s u r + cu r - 7o|x| a ~ 1 u r = Cu p in V'(B N (0, r)). 
We have, by (5.20), 

(-A) s u r - jalx^ur > (( 7o - ~i a )r a ~ l - c)u r + Cu p > Cu p in V'(B N (0, r)). 
Therefore by Lemma 3.2, there exists a constant C r > such that 

(5.22) u r > C r \x\-^ +a inB N (0,r/2). 
For every fceN, take v k G J% s (B N (0,r)) as the solution to 

B s v k + cv k - ( 7o - l/k)\x\ a - l v k = Cmin«, k) in 5^(0, r). 
Since u r satisfies (5.21), it follows (see Remark 2.9) that 

(5.23) u r >v k >0 in B N (0,r) 
by Lemma 2.3 and Lemma 2.2. 

Put V k =U{v k ) we have that for any ^ G H^ T (B N (0,r);t a ) 
(5.24) 



k- 1 t a VV k -V^dz = (jo-l/k) \x\ a ~ l v k ^dx- c i v k Vdx 

Jr^ +1 JB N (0,r) JB N (Q,r) 

+C I min(u p ,k)^dx. 

JB N (0,r) 

Thanks to (5.22), we can choose r' G (0,r/2) (small) such that 
(5.25) -c + CvP' 1 > 1 in B N (0, r'). 

For such a fixed r', take ip G C™(B N (py)). For e > 0, set l/ fc £ = + e and put 
V> := 2^1. We have 

|V^M| 2 = |V£VVf + W fe • v(ifv 2 ) 
and also V^ 2 G H^ T (B N (0,r);t a ). This together with (5.24) yields 

k- 1 ! t a \vn(<p)\ 2 dxdt > (70-1/fe) / isr -1 — ^— 

Jr^ +1 JBN(oy) v k + e 



B N (0,r') 

ip 2 dx 



B N (0,r') v k + £ 
f min(? 
B N (0,r') v k + ^ 



.//>• 
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Thus 



1 / t a \VU{<p)\ 2 dxdt > (70-l/fc) [ Ix^-^—^dx-c f <p 2 dx 

Jr i ?+ 1 JB N (o,r>) v k + e jB N (0,r') 



r (0,r') v k "I" £ JB N (0,r>) 

+C [ min( ^ fc V dx. 

JB N (0,r') v k+£ 

By Fatou's lemma, when e — > 0, we have 



kJ 1 / t a \Vn{f)\ 2 dxdt > (70-l/fc)/ \x\ a - x ip 2 dx- c \ f 2 dx 

Jr1 +1 JB N (oy) JB N (0,r') 



B N (0,r') JB N (0,r') 
mm(u?,k) 2 



+C [ ""^"'"V ^. 

JB N (0,r>) v k 



By (5.23) we obtain 



1 / t a |V?%)| 2 dxcit > (70-l/fc)/ \x\ a - l <p 2 dx - c [ (p 2 dx 

JR N+1 JB N (0,r') JB N (0,r>) 



+C I min( ^' fc) ^dx. 

JB N (0,r') u r 

It follows again from Fatou's lemma that 

k; 1 / t a |V^(v?)| 2 dxdt > 70 / |x| a "V 2 dx+ / (-c + Cur 1 )^ 2 ^- 

JR^ +1 JB N (0,r') JB N (0,r') 

From the choice of r' in (5.25), we get immediately, for any ip G C^°(S Ar (0, r')), 
k; 1 / t a \VH(ip)\ 2 dxdt-j [ \x\ a - 1 i P 2 dx> [ p 2 dx. 

JR^ +1 JB N (0,r>) JB N (0,r') 

By scaling we have (5.19) which was our objective. 

□ 
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